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We consider a version of QCD dipole cascading corresponding to a finite number n of discrete AY 
steps of branching in rapidity. The disctretization scheme preserving the holomorphic factorizability 
and scale-invariance in position space of the dipole splitting function, we derive an exact recurrence 
formula from step to step which plays the role of a "discrete-time" Balitsky-Kovchegov equation. 
The BK solutions are recovered in the limit n = oo and AY = 0. 
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(N 1 I- INTRODUCTION 

X) _ 

D ' As was first suggested in [lj, the distribution of partons inside a hadron, in the limit of high-energy, and fixed 
(large) Q 2 and large N c , can be approximated by a cascade of colourless dipoles. Starting from this observation, and 
£NJ ' using the evolution equation describing the development of the cascade in rapidity Y, it was possible to obtain -in 
the leading logarithmic approximation- an evolution equation for the scattering amplitude of the dipole cascade on 
an "uncorrelated" target satisfying the unitarity constraint Q • The dipole cascading formalism gives a derivation of 
QCD evolution equations obtained from the perturbative expansion in the leading logs approximation in energy. Both 
the linear regime corresponding to the Balitsky Fadin Kuraev and Lipatov (BFKL) evolution equation and its 
non- linear extension, the Balitsky-Kovchegov (BK) equation @, Bj, find a convenient description in terms of dipoles. 
A general characteristics of the solutions of the BK equation was recently discussed in J5j. It was shown that an 
approximate geometrical scaling is the generic asymptotic energy property of the system, related to mathematical 
solutions in terms of traveling waves, and a general method of finding these solutions was developped. 

The dipole cascading formulation of Refs.0, Q corresponds to a classical branching process where dipoles split 
with a probability distribution given by the BFKL kernel [!| expressed in transverse position space taking the role 
■ of (2-dimensional) space, the rapidity variable having the role of time. As such it belongs to the large family of 
Oh] random branching processes which is largely studied in statistical mechanics and in mathematics (see |a, la, 13 la l2||)- 
While the original QCD dipole formulation of [J, 0] is based on random branching, we remark that well-known classes 
| of branching processes are considered with discrete steps in time ||, with interesting physical and mathematical 
properties, and applications. We want to address here the question of a similar discretization of dipole cascading. 
The attractive feature of the QCD dipole approach is that -being formulated in the framework of the QCD per- 
. , turbation theory- it is a well-defined stochastic fragmentation system. In particular the cascade vertices are uniquely 
given by the theory. Any departure from this scheme runs into a problem of serious ambiguities even if one imposes 
J~j ] the condition that in the some limit one should recover the known perturbative results. 

In the present paper we study one class of dicrete-in-rapidity dipole cascades whose vertices are scale-invariant 
in transverse space. In particular we retain the important feature of holomorphic separability which is present in 
the leading-logarithmic approximation. We show that such a theory naturally leads to discretization of the dipole 
cascade and that in the limit of a large number of rapidity steps one recovers the results of the leading-logarithmic 
approximation. One may hope that this exercise will help to understand the structure of gluon cascading beyond the 
lowest order of the perturbation theory. 

The plan of our study is the following. In the next section 2, we derive a discretized version of the dipole cascading 
preserving scale-invariance and the holomorphic separability of the original BFKL kernel. In section 3 we derive the 
corresponding non-linear master equation, which plays de role of a "discrete time" Balitsky-Kovchegov (discrete-time 
BK) equation. In the next section 4, we use the method of traveling wave solutions in order to explore the solutions 
of the discrete-time BK equation when the number of steps goes to infinity. In particular we prove the convergence to 
the BK solutions when the rapidity step size goes in the same proportion to zero. In section 5, we draw some possible 
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FIG. 1: Discretized branching process of QCD dipoles. The figure shows the rapidity steps for Y evoiution of one dipole xo, xifor 
successive splittings at transverse space points X2 and then 2:3,3:4, etc... 



interesting outcome of our discretisation procedure. 



II. THE MASTER EQUATION 



Starting from the the Balitsky-Kovchegov equation we show how the modification of the vertex leads to the dis- 
cretization of the dipole cascade. 
The BK equation reads: 



dS f 

— (x i,Y) = / d 2 x 2 



K(x ,x 1 ;x 2 ) {S(x 2,Y)S(xi2,Y) - S(x i,Y)} , 



(1) 



where S(xqi,Y) is the S-matrix element for the dipole-target amplitude, xqi — xo — xi being the dipole transverse 
size. The BFKL kernel 



K.(xq, xi;x 2 ) d 2 x 2 = 



27T Xr\nX I 



d 2 x 2 



(2) 



t 02 a '12 



has the classical interpretation fH of a splitting probability in transverse space by unit of rapidity. Note that this 
kernel verifies simultaneously scale invariance, the X02 <-> X12 symmetry and holomorphic separability. To see this 
last property, it is worth using the 2-dimensional components of scaled variables £^ = Xi/xqi and transforming into 
complex variables, namely 



(i) 



- ^ 



(2) 



(1) 



(2) 



(3) 



Holomorphic separability means the factorization property of the kernel J5J into analytic (function of z) and antian- 
alytic (function of z) parts. It plays an important role in the conformal symetry properties of the BFKL amplitudes 

m 

Let us consider a modified kernel preserving this property: 



K,(xq,Xi;x 2 ) d 2 x 2 



a.N B 



2lT 2 \2:o2 :I '12 



d 2 X2 _ ct s N c 



'01 



2tt 2 



(zz(l-z)(l-E)) 



a—i dzdz 

"IT' 



and call 



N = / d 2 x 2 jC(x Q ,xi]x 2 ) , 



(4) 



(5) 



3 



where TV" is finite for < a < 1/2. 

We can now safely work separately on the two terms of Q and write 

JJ~^ = jy d 2 x 2 JC d (xo,x 1 ;x 2 ) S(x Q2 ,Y)S(x 12 ,Y)^ - S(x i,Y) , 

where Kd = K./Af is now a properly normalized probability distribution. 
Using the known mathematical identity 1 

— [ dzd - z Z ^- Z (^) (i - *r-v (i . 4^) = mjm ni-A-B) 

2m J y ' K ' r(i - a) r(i — b) t(a + b) 

we obtain the expression for the normalization 



and the following expression for the probability distribution 

d 2 x 2 v f 

r 2 (l-a)r(2a) 



7T \ X r\i-> X 1 



(6) 



(7) 



_ a s /Y c r 2 (a) r(l - 2a) 

^ " "17" x r 2 (l-a)r(2a) ' (8j 



. o if < \ 1_a d 2 x 2 fr 2 (a) r(l-2a)". 



Choosing a finite interval 

TVAF = An = 1 (10) 
and considering neN,ti>l allows one to transform Eq. © into a finite difference equation 

5„+i(.T i) - S n (xoi) = J d 2 x 2 JCd(xo,x 1 ;x 2 ){S n (x 2)S n (x 1 2) - S n (x i)} , (11) 

leading to (the distribution ICd, see ||SJ|, being normalized to 1) 

S n+1 (x i) = / /Cd(x ,a;i;x2) d 2 x 2 x S„(iro2) 5„(xi 2 ) . (12) 



As clear enough from its formal structure, Eci. (|12fl is a "discrete time" version of the Balitsky-Kovchegov equation 
i ffl l. It remains to be proven that it leads back to the BK equation when going to its continuous limit, which will be 
defined in the next section. 

Indeed, Ea.lfT2)l has the typical structure for S-matrix elements defined for a branching process (or tree structure) 
with discrete steps of time evolution, see Fig^ At each step n — ► n + 1 a dipole of size xqi splits into two dipoles 
of sizes xo 2 ,xi 2 at the point x 2 with a 2-dimensional probability distribution JCd(xo, x\\ x 2 ). In this description, the 
formula l|l(J|) determines the length of the "rapidity veto" AY = 1/N . The dependence of AY on a is shown in Figure 
2. 



III. THE RECURRENCE STRUCTURE IN MOMENTUM SPACE 

An even simpler form of equation (|12fl can be obtained in momentum space, when one considers solutions indepen- 
dent of the impact-parmeter. It reveals even better the recurrence structure of the discrete-time BK equation. Using 
a 2-dimensional Fourier transform, one defines 

= J^(4i) a e^S n (x al ) J dx j o(kx) Sn (x) 

" [ '~ J (sgja e ik.x 01 J dx J (kx) X 2 "- 1 ' K } 



1 There is one condition for to be valid, namely A — A, B — B S Z, which also ensures that formula is symetric in the interchange 
A -* A,B -> B. 
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with 

/ dx J (kx) x 2 "- 1 = kr 2a -^^ . (14) 

J r (!- a ) 

working in the context of an impact parameter independent framework. Note that S„ (k), can be obtained from a 
convolution of the S-matrix element in momentum space by the Fourier transform of the the weig ht (x^) in the 
integrand of 

It is straightforward to infer from (|12|) and the following relation 

CIO*) = (^TV-) 2 I ^ ^(*o 2 ) Q {^T e*-(-^ S n {x m ) S nM , (15) 

V \ J / J x 12 X 12 

where the substitution of the integration variable x 2 — > x%2, with Jacobian unity, has been performed. 
Then, the discrete-time BK equation (|12ll can be rewritten in a particular simple way as 

sffi(k) = {s<°>(k)} 2 , (16) 

which clearly expresses the nature of the non-linear recurrence relation from step to step, this recurrence structure is 
also clear when expressed in terms of transition matrix elements with the same conventional notation as l|13|l 

f^(k) = 2f(«Hk)-{f^(k)} 2 . (17) 



IV. HIGH ENERGY LIMITS 
A. Linear regime 

Let us examine the properties of the master equation l|12|) (also (|1€>|) . Considering first the linearized form of Eq. 
(|12fl near S ~ 1 we have 

Xi+iix^) = J d 2 x 2 IC d {x ,x 1 ;x 2 ) {%i(xl 2 ) + %i( x i2)} > ( 18 ) 
where T = 1 — S is the transition matrix element. Using the complex scaled variables J3J one writes 

Zi+i{t) = 2 J dzdz fC d (z,z) T n (tzz) , (19) 

where 

( r 2 (i- Q )r(2 Q ) \ ! 
^'^ U 2 ^) r(i- 2o ) J W-^a-^)} ■ CM) 

and the factor 2 comes from the z 1 — z symmetry of K-d- 
We now introduce the Mellin transform representation 

Tn{t)= Sc^t 11 fn{l) ' (21) 

and easily derive the recursive relation 

f n+1 ( 7 )=T„( 7 ) e* (a) W (22) 

where, using Q, 

X<«»(7) . log (2 r r 'Y r ?' r(1 -^' 1 - 2 r r ?2?1 ) ' < 23 > 
\ 1 (a) 1 (1 — a — 7) 1 (1 — 2o) 1 (2a + 7) / 



Hence, the solution of the linearized equation is simply 

t n ( 7 ) = To (7) exp {n X (a Hl)} = Ml) exp [MY x [a) (l)} , (24) 

where Y is the total rapidity interval and M is given in (JSJ. 

Let us now consider the limit a — > 0. A straightforward algebra gives 



N x {a) = { ~ + 0(a) \ x {a [2*(1) - *( 7 ) - ¥(l- 7 )] + ^(a 2 )} 



^( 2 tf(l)_tf( 7 )-*(l- 7 )) 



where one recognizes in the brackets the BFKL kernel function x( 7 ). Hence, from i|24|) 



■0(a) - 
(25) 



T n ( 7 )=To( 7 )exp{no[(2*(l)-*(7)-^(l - - 1] ) =%[') <-x P <j^x(-> > (-><!> 
in which we recover the solution of the BFKL equation. 



B. Non-linear regime 

Let us now consider, in the spirit of [5j for the Balitskii-Kovchegov equation, the asymptotic solutions of the master 
equations (|12|l or in the continuous limit, i.e. for a large number n of steps of the cascade. Asymptotic solutions 
we are discussing now may be valid only at really extreme energies. Nevertheless we feel that they may be of some 
interest. Indeed, as we shall see now by inspecting the general properties of the full non-linear equations for arbitrary 
a, the master equations (|12|l (or equivalently {TSJ) give the same solutions at large n than the Balitsky-Kovchcgov 
equation in the limit a — > or, in more mathematical terms, stay in the same universality class. We will also extend 
our study to another type of continuous limit, namely when n is large but a is kept fixed. 

Let us start from the expressions (|21|l and (|24|l giving the solution of the linearised problem. One has 

% ® = JJ^ t 1 % (7) exp [MY x [a) (7) } , (27) 

where x has been obtained in l|23() . We will now follow mathematical arguments used for the solution of discrete 
non-linear equations appearing in a statistical mechanic context 0. They are similar than those used in 5] for the 
asymptotic solutions of the Balitsky-Kovechegov equations and are part of more general mathematical results on 
non-linear equations, such as for the Fisher and Kolmogorov, Petrovsky, Piscounov equations 0,0 of "pulled front" 
type 0. 

Similarly to the arguments developped in |5(, the solution of the non-linear equation equations (|12|l or are 
traveling waves whose expressions can be obtained starting from the linear equation l|27|) . 
Let us reinterpret l|27|l as a linear superposition of waves: 

T n (t)= f c §^?oh) exp{- 7 (L WF + H+^X (a) (7)} , (28) 

where L = logl/t and L WF = L — vn is looked for as the scaling variable determining the region moving with the 
traveling wave front, x^ (7) defines the dispersion relation of the linearized equation. In particular, each partial wave 
of wave number 7 has a phase velocity 

M7) = (29) 

7 

whose expression is found by imposing that the exponential factor in Ea. (|28J) be independent of n for v = v v (j). By 
contrast, the group velocity is defined by the saddle point 7 * of the exponential phase factor 



The key point of the mathematical derivation of the asymptotic solution of the non-linear equation is that, for 
appropriate initial conditions 0.0. the critical regime at 7 = 7 * is selected by the non-linear damping. 




The velocity v of the front is denned by 



X W( 7 *) . X W( 7 ) 
i; = = mm 

7* 77 

where the value of 7* is determined from the equation 

dx (a) x (a)( 7) 



dr) 
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(31) 



(32) 



7* 



Indeed, in this case, the group velocity is identical to the minimum of the phase velocity 2 . Now, since 

n = Y/AY = YAf (33) 

we can write v*n = vyY , with 

v Y = v* TV. (34) 
Reporting these results in Ea. (|28|) gives the dominant term in the asymptotic expression for the form of the front: 

T(y,t)~To(7*) exp{-7*loglA + /vy Q )( 7 *)y} . (35) 

Note that two more "universal prefactors" (ie. independent of the initial conditions and of the precise form of the 
non-linear damping terms) can be obtained 0,0] m the asymptotic expansion of the solution. 

From H35J) . it is clear that the same result than for the Balitsky-Kovchegov equation will be obtained in the limit 
a — > 0, due to the equivalence of the linear kernels and the universality properties due to the the non-linearities. This 
achieves the proof. 

Let us consider our high energy limits for n — > 00 with a ^ 0, 1/2 kept fixed. Some features of the asymptotic 
solutions are displayed in Fig[21 As shown in the left part of the figure, the rapidity step Ay is displayed as a function 
of a for a given value of a s = .1. There is a maximum of the reachable rapidity step for a ~ .3 whose value depends 
on a s . At this stage this remains an intringuing feature of our scheme of discretization. The right hand part of 
the figure gives the effective slope (or "intercept") in rapidity, which, starting from the BK value stays constant till 
approximately a = .25 and then grows. Note that the energy momentum constraint on the number of steps at given 
rapidity range cannot be taken into account in this continuous limit, as will be discussed in the next section. 



2 The relation v g = v v , in analogy with wave physics, has been written in Ref. 
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V. CONCLUSION AND OUTLOOK 



To summarize the main results of our paper, we have considered a "discrete time" version of QCD dipole brancing 
process characterizing the approximation of the perturbative expansion in the leading logs of the energy. This 
is obtained by a modification of the splitting probability for finite steps in rapidity preserving the holomorphic 
factorizability of the initial BFKL kernel. 

Our main result is the derivation of an exact counterpart of the BK saturation equation in terms of recurrence 
formulae of quite simple form, expressed both in position space l|12|) or in a (modified) fourier-transformed space I jlfijl • 

We checked that we recover the BFKL and BK equations in the appropriate AY — > continuous limit while the 
high energy asymptotics for AY ^ have been investigated. 

The principal outlook of our study concerns the problem of subasymptotics. Indeed, the master equations Ijl2|l 
and Ijl6(l allow for an iterative solution of the non-linear evolution problem. Inserting any initial condition 7o in 1|12|) 
or in (modified) Fourier transform for (|16|l , it is possible to generate the full solution, at each step of evolution. In 
particular, as explained in Section 2, for a ^ 0, 1/2 the steps of the cascade are separated by a finite distance AY 
in rapidity. Therefore, for a given total energy, the number of steps in the cascade are limited by n max « Y/ AY 1 
where Y is the total available rapidity. Hence these equations give a convenient way to examine the effect of a limited 
number of steps, and thus of dipoles, in the cascading process. Indeed, the gluon cascade description of Q, can 
only be justified in the leading logarithmic approximation. The investigations of higher order corrections shows pl| 
that they tend to limit the number of emitted gluons 3 . One possibility to take this effect into account is to forbid 
the emission of gluons which are too close in rapidity 0, Q] . Since the emitted gluons are at the source of dipole 
splitting [l| , it is expected that the formation of dipoles is similarly limited in rapidity. Such a "rapidity veto" can 
-in turn- be approximated by a discretized cascade where gluon (or dipole) emissions is separated by a finite distance 
in rapidity. We thus feel that our equations (|12I16|| ) could give an efficient way to investigate this problem. It also 
can be useful to develop the parallel with statistical mechanics properties which appeared to be recently fruitful @ . 

When completing the writing of the present work, the paper [l5j has appeared, claiming that the solution of an 
evolution with discrete rapidity intervals may lead to a chaotic behaviour. It stems from a toy model of the BK 
equation in the form of a logistic map. We note that our equation (fTTfl has the structure of a logistic map if one can 
neglect the shift (a) — > (2a) in the weight factor. However the "Malthusian parameter" would be then 2 instead of 
3.77 (see the discussion in Ref.[l5j) and thus is smaller than the value leading to chaos. Since our equation l|17(l is a 
consistent discretization of the BK equation, the investigation of its convergence (or not) to the mean fied problem 
deserves to be studied. 
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This is partly due to the simple effect of energy conservation which is ignored in the LL approximation. 



